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§1  Introduction 

In  this  report  we  present  a  block  diagonalization  theorem  which  is  designed  to  study  the 
stability  and  bifurcation  of  rotating  systems,  or  more  generally,  of  relative  equilibria.  The  context 
of  the  discussion  is  the  energy-momentum  method  for  mechanical  systems  with  symmetry.  Simo, 
Posbergh  and  Marsden  [1989]  and  Lewis  and  Simo  [1989]  discovered  crucial  special  cases  of  the 
block  diagonalization  theorem  for  uniformly  restating  systems,  including  general  nonlinear  elasticity 
and  geometrically  exact  rods.  Our  purpose  is  to  abstract  these  examples  and  prove  a  general 
^  geometric  theorem.  We  expect  these  general  results  will  be  important  for  rotating  gravitational 
I _ fluid  masses  as  well. 

For  rotating  systems  the  result  says  that  a  splitting  of  coordinates  can  be  explicitly  found  on 
a  linearized  level  which  represent  the  rotational  and  internal  vibrational  modes.  In  these 
coordinates,  the  second  variation  of  an  augmented  Hamiltonian  is  block  diagonal.  Of  course 
coordinates  can  always  be  found  in  principle  to  do  this,  but  we  are  able  to  do  it  explicitly  enough  to 
give  useful  stability  and,  we  also  believe,  bifurcation  criteria. 

On  the  other  hand,  the  symplectic  form  does  not  block  diagonalize,  indicating  that  the 
rotational  and  internal  modes  are  in  fact  dynamically  coupled.  However,  for  purposes  of  the 
stability  calculation,  block  diagonalization  of  the  augmented  energy  is  what  is  important.  The  off 
diagonal  terms  in  the  symplectic  form  (sometimes  called  Coriolis  coupling  terms)  are,  however, 
sufficiently  simple  that  they  should  be  useful  for  studying  the  dynamic  interaction  of  the  rotational 
and  internal  vibrational  modes. 

For  rotating  pseudo-rigid  bodies,  Lewis  and  Simo  [1989]  noticed  that  the  computation  of 
the  definiteness  of  the  second  variation  is  considerably  simplified  by  our  result  -  in  this  case  the 
simplification  saves  considerable  computation  time.  In  fact  the  symbolic  and  numerical 
manipulation  required  one  to  test  a  full  14  x  14  matrix  for  definiteness;  block  diagonalization 
reduces  this  to  testing  a  6  x  6  matrix  for  nonisotropic  bodies  and  to  a  3x3  matrix  for  the  isotropic 
case. 


§2  The  Energy-Momentum  Method 

We  begin  our  work  in  the  context  of  standard  mechanical  systems  with  symmetry  before 
any  reductions  have  taken  place.  In  other  words,  we  begin  with  a  symplectic  manifold  (P,  Q) 
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rather  than  a  Poisson  manifold.  In  fact,  shonly  we  shall  specialize  to  the  case  of  P  =  T*Q  and  a 
Hamiltonian  of  the  form  kinetic  plus  potential. 

Let  G  be  a  Lie  group  acting  symplectically  on  P  with  an  equivariant  momentum  mapping 

J:P^8*  (1) 

(see  Abraham  and  Marsden  [1978],  Marsden  [1981]  or  Marsden  et  al.  [1982]  for  the  standard 
definitions  and  results  used  here). 

Let  H  :  P  R  be  a  given  G-invariant  Hamiltonian.  A  point  in  P  is  called  a  relative 
equilibrium  if  there  is  a  ^  e  8,  the  Lie  algebra  of  G,  such  that  for  all  t  e  R, 

z(t)  =  exp(t4)  Zg,  (2) 

where  z(t)  is  the  dynamical  orbit  of  the  Hamiltonian  vector  field  of  H,  with  z(0)  =  z^. 

The  energy-momentum  method  rests  on  the  following  result. 

2.1  Relative  Equilibrium  Theorem  A  point  z^  is  a  relative  equilibrium  iff  there  is  a  ^ 
€  8  such  that  z^  is  a  critical  point  o/  H^ :  P  ->  R,  where 

H^(z)  =  H(z)-<J(z)-p,,^>  (3) 


and  n^  =  J(Zg). 


In  (3),  the  Lie  algebra  element  ^  €  8  may  be  regarded  as  a  Lagrange  multiplier.  Since  J 
is  conserved  by  the  flow  of  Xjj,  the  set  J  -  =  0  is  preserved,  so  one  may  regard  it  as  a  (non- 

holonomic)  constraint  set.  It  also  follows  that  ^  e  8n_.  the  isotropy  algebra  of  p  (with  respect 

He 

to  the  coadjoint  action).  Thus, 


6H^(z^)  =  0  (4) 

may  be  regarded  as  a  (constrained)varjari<?na/  principle  for  relative  equilibria. 

The  relative  equilibrium  theorem  is  readily  verified.  Of  course  it  has  a  long  history,  going 
back  to  Lagrange  and  Poincard  for  rotating  systems.  Like  many  basic  results,  it  has  been 
rediscovered  in  a  number  of  contexts  by  various  authors.  Early  references  in  our  context  are 
Arnold  [1966],  Smale  [1970]  and  Marsden  and  Weinstein  [1974]  As  we  shall  state  below,  the 
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relative  equilibrium  theorem  sometimes  specializes  to  the  principle  of  symmetric  criticality 
(Palais  [1979]). 

The  energy-momentum  method  proceeds  as  follows  (see  Holm  et  al.  [1985]  for  the 
meaning  of  formal  stability  and  related  references). 


Energy-Momentum  Method 


To  test  a  relative  equilibrium  €  P  for  formal  stability: 

1  Choose  q  such  that  5H^(2g)  =  0  T 

2  Choose  a  linear  subspace  S  c  T^^P  such  that 
i  S  c  ker  dj(z J  and 

ii  S  complements  T  (G  •  z  )  in  kerTJ(z  ),  where  G  c  G  is  the  isotropy  — 

He  ^  "  He  -  -■ 


The  energy-momentum  method  "covers"  the  energy-Casimir  method  (Holm  et  al.  [1985]) 
in  the  sense  that  if  the  latter  applies  and  gives  formal  stability,  so  does  the  former.  One  difficulty 
with  the  energy-Casimir  method  is  that  on  the  reduced  space  P/G,  there  may  not  be  enough 
Casimirs  to  make  the  method  effective,  even  to  get  the  analogue  5(H  +  C)  =  0  of  (4).  This 
difficulty  is  genuine  for  the  case  of  geometrically  exact  rods,  for  instance.  See  Simo,  Posbergh 
and  Marsden  [1989]  for  further  details. 

The  fact  that  8^^(Zg)  drops  to  the  reduced  space  follows  from  the  next  lemma. 


2.2  Gauge  Invariance  Lemma 


62H^(z,)(Tip(z^),5z)  =  0 


(5) 


for  all  5z  e  ker  TJ(Zg)  and  i]  e  g  ,  where  Tjp  denotes  the  infinitesimal  generator  of  the  group 
action  on  P. 
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This  follows  readily  from  invariance  of  H  and  equivariance  of  J.  One  can  view  (5)  as  a 
block  diagonalization  result  on  the  unconstrained  tangent  space  T  P,  but  it  does  not  yield  block 

diagonalization  within  the  constrained  subspace  S  in  the  energy-momentum  method.  It  is  the 

latter  that  we  are  concerned  with. 

One  can  identify  any  choice  of  S  with  the  tangent  space  to  the  reduced  space 

at  [Zg]  (assuming,  as  we  shall,  that  |ig  is  a  regular  and  generic  value;  c.f.  Weinstein  [1984]). 
However,  it  is  easier  to  do  our  analysis  direcdy  on  T^T*Q  rather  than  on  the  quotient  space.  This 

is  the  usual  situation  found  in  constrained  optimization  problems.  However,  dropping  the 
calculations  to  the  quotient  space  at  the  appropriate  point  will  play  a  useful  role. 


§3  Simple  Mechanical  Systems 

We  now  specialize  to  the  systems  we  will  be  studying.  Let  Q  be  a  configuration  manifold 
and  P  =  T*Q  with  its  canonical  symplectic  structure  and  cotangent  coordinates  (q',  Pj)  in  the 
finite  dimensional  case.  {Whenever  we  use  coordinates,  we  assume  Q  is  finite  dimensional, 
although  the  results  are  mt  restricted  to  this  case.)  Coordinates  on  the  velocity  phase  space  TQ 
are  similarly  denoted  (q‘,  q*). 

Let  g  denote  a  Riemannian  metric  on  Q;  in  coordinates  we  write  the  components  of  g  as 
gjj  as  usual,  and  we  write  g‘J  for  the  inverse  tensor.  Let  K  :  TQ  — >  R  denote  the  corresponding 

kinetic  energy,  i.e., 

K(q,q)  =  5-gy(q)q'^  .  (D 


Let  V  ;  Q  R  be  a  given  potential. 

Assume  G  acts  on  Q  (by  a  left  action)  and  hence  on  T*Q  by  the  cotangent  lift,  so  the 
equivariant  momentum  map  is  given  by 

<J,^>(aq)  =  (aq,^(q)).  (2) 

In  coordinates,  we  define  the  action  coefficients  A^(q)  by  writing 
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[^(q)]'  =  A'(q) 


(3) 


where  a,  b,  c, ...  denote  coordinate  indices  for  the  Lie  algebra  9.  Thus  (2)  becomes 

Vq.P)  =  PiAi(q).  (4) 

We  assume  that  G  acts  on  Q  by  isometries  and  that  the  potential  V  is  G-invariant.  For 
elasticity,  for  instance,  this  is  the  requirement  of  material  frame  indifference.  Note  that  (3)  of  §2 
reads 

H^(q,  P)  =  ^  g'j  Pi  Pj  +  V(q)  -  pj  A^(q)  (5) 


Define  the  moment  of  inertia  tensor  I  for  the  system  locked  at  qe  Q  by 

lab(q)  =  gij(q)  Ai(q)  A^(q)  (6) 

(alternatively,  in  terms  of  the  q-dependent  inner  product  (^,  Tj)  :=  {^(q),  'nQ(q))  on  g,  we  have 
(^,  T|)  =  I,b(q)^*T|*’).  and  define  the  augmented  potential  by 

V^Cq)  =  V(q)-ll.t(q)4*^‘'.  (7) 

One  can  readily  verify  the  following  (see  Abraham  and  Marsden  [1978]  and  Palais  [1979])  by 
writing  out  the  conditions  5H^  =  0  in  2.1 .  A  more  elegant  argument  is,  however,  given  below. 

3.1  Principle  of  Symmetric  Criticality  A  point  =  (q\  pj)  is  a  relative  equilibrium 
if  and  only  if  there  is  a  9^1^  such  that 

I  Pi  =  gjj  Aj  (i.e.,  Pg  is  the  Legendre  transform  of  ^Q(q^))  (8a) 

and 

g 

ii  q^  is  a  critical  point  of  — Vt  I  =0.)  (8b) 

^  6q  ^  I  q  =  qe 

This  is  useful  for  carrying  out  the  computations  that  follow.  We  also  observe  that  is 
G^-invariant,  and  so  induces  a  function  on  Q/G^. 

Define  the  one-form  A^  on  Q  by 
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A?(q)  =  gy(q)  AJ(q)  (9) 

or  abstractly,  A^q)  =  [^(q)]^  ,  where  **  denotes  the  index  lowering  operation  with  respect  to  the 
metric  gjj.  In  other  words,  A^(q)  is  the  Legendre  transform  of  ^gfq).  We  remark  that  A  may 
be  viewed  as  a  G-connection  for  the  bundle  Q  Q/G  and  that  this  connection  plays  an  important 
role  in  Berry's  phase;  cf.  Marsden,  Montgomery  and  Ratiu  1 1988].  Now  notice  that  at  equilibrium, 
(8a) says 


Pe  =  A^(qe)  ■ 


(10) 


Also  note  that 


H^(q,p)  =  K^(q,  p)  +  V^(q,  p)  +  (p^,  0  (11) 

where  K^(q,  p)  =  j  ||  p  -  A^(q)  |p  ,  and  is  given  by  (7).  By  (10),  has  a  critical  point  at 
Zg.  Thus,  (8b)  is  a  direct  consequence  of  the  relative  equilibrium  theorem  and  (11). 

In  the  energy-momentum  method  we  shall  use  a  special  choice  of  S,  namely 

=  { ''ze  ^  I  TtIq  •  v^  is  g-orthogonal  to  T(G^^  •  q^)  and  v^  e  ker[TJ(Zg)] ) .  (12) 

Letting  coordinates  on  TT*Q  be  denoted 

(q'.  Pi.  5q'.  Spj), 


(12)  reads,  with  the  help  of  (8a), 

.5  =  { (q',  Pi,  8q',  Spj)  I  gjj(5q)' AJ  x“  =  0  for  all  X  ^  3^^  and 
(5p)i  A*  -I-  gjj  A^  ®  } 


(12’) 


§4  Rigid  Variations 

One  version  of  the  cotangent  bundle  reduction  theorem  (see  Abraham  and  Marsden  [1978] 
and  Kummer  [1981],  Montgomery  [1986]  and  references  therein)  states  that  the  reduced  space 
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(T*Q)^  is  a  symplectic  bundle  over  T*  (Q/G)  with  fiber  the  coadjoint  orbit  through  |ig.  Thus 
there  is  an  isomorphism 

9/9e  X  Tf^(T*(Q/G))  =  g/9^x  (Uj^^  xV*^) 

where  is  a  model  space  for  Q/G.  For  G  =  S0(3),  Ujjsfp  models  the  configuration  space 
for  the  internal  modes,  while  models  the  phase  space  for  rigid  modes.  Our  goal  is 

to  realize  this  decomposition  explicitly,  in  such  a  way  that  8^H^(z^)  block  diagonalizes.  The 
bundle  (T*Q)^ ->  T*(Q/G)  with  fiber  also  has  a  natural  connection  (Montgomery  [1986]) 

and  our  decomposition  should  be  related"  in  some  way  to  the  horizontal-vertical  split  for  this 
connection.  However  we  proceed  directly  here;  see  also  the  comments  in  §5  below. 

We  will  define  two  subspaces  i^jQ  and  of  S  and  further  subspaces  and 


“  ‘%IG  ®  ‘^INT  =  *^10  ®  (^INT  ® 

relative  to  which  5^H^(Zg)  will  be  block  diagonal.  As  above,  the  first  component  =  8/9^^ 
of  S  is  isomorphic  to  the  tangent  space  to  the  coadjoint  orbit  through  p.^.  As  we  shall  see,  this 
component  will  also  carry  the  coadjoint  orbit  symplectic  structure.  This  first  component  is  defined 
in  terms  of  rigid  variations  as  follows:  Let 

9q  =  {TlQ(q)e  TQ1ti€  gandqe  Q)  (2) 


and  let  TggcTTQ  be  its  tangent  bundle. 


4.1  Definition  Let  Vj^jQ  =  sfTgg)  where  s  ;  T^Q  — >  T^Q  is  the  canonical  involution. 
Alternatively,  Vj^jq  consists  of  double  tangents  of  curves  denoted  by  Aq  {identified  with 
velocity  variations  of  superposed  rigid  body  motions  in  the  case  of  S0(3)) 


where  T[(t)  is  a  curve  in  g  with  'n{0)='n  and  q(t)  is  a  curve  in  Q.  (The  canonical  involution 
in  effect  swaps  the  order  of  differentiation.) 
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In  coordinates,  if  we  write  elements  of  Vj^jq  as 

(q*.  q‘,  Aq^  Aq‘).  (3a) 

then  we  find  that 

dA' 

Aq'  =  and  Aq‘  =  — q’^  t)®  +  A^  (3b) 

dq*' 

Now  let  FL :  TQ  -» T*Q  be  the  Legendre  transform  given  by 

Pi  =  gijO^  W 

and  let  TFL  :  TTQ  — >  TT*Q  be  its  tangent  map.  Set 

‘^RIG  “  ■  ^RIG 

where  S  is  defined  by  (12)  of  §3  .  If  we  let  a;f  denote  the  (q-dependent)  orthogonal 
complement  of  in  the  metric  then  one  finds  that  is  parametrized  by  elements  T)  e 
a^  as  follows:  we  write  elements  of  as 


(q',  pj,  Aq‘,  Apj) 

(6a) 

where 

Aq*  =  A;ti® 

(6b) 

and 

dAl^ 

APi  =  --^^^‘Pk+gijAlC" 

(6c) 

where  Tj  €  ai.  and  where  ^  e  a;  the  condition  that  (6a)  belongs  to  ker(T  J)  is  equivalent  to  the 
relation 

;■  =  I">(ad;  (7) 

i-e..  ?  =ad;n,,  SO  ^  is  determined  by  T).  One  checks  that  t)  e  a^J^  as  well. 
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§5  The  Internal  Vibration  Space 

Now  we  define  a  complement  to  ^iq  in  S.  We  will  do  this  by  a  constructive  procedure 
that  can  be  effectively  carried  out  in  examples.  As  we  have  mentioned,  this  complement  appears  to 
be  not  the  same  as,  but  related  to  the  complement  to  the  vertical  space  relative  to  a  natural 
connection  on  the  coadjoint  orbit  bundle  (T*Q)^  — >  T*(Q/G).  In  this  regard  we  note  that  the 

metric  naturally  induced  on  Q/G  is  Wilson’s  G-matrix  (see  Wilson,  Decius  and  Cross  [1955]). 
Our  decomposition  appears  to  be  finer  than  the  one  proposed  by  Guichardet  [1984]  and  discussed 
by  Iwai  [1988].  Notice  that  we  have  connections  on  all  levels  of  this  tower  of  bundles 

(T*Q)^  T*(Q/G) 

where  >  (T*Q)|j^  is  regarded  as  a  G^  bundle  and  (T*Q)^  — >  T*(Q/G)  is  regarded  as  an 

bundle,  where  is  the  coadjoint  orbit  through  fx. 

The  Guichardet-Iwai  results  appear  to  be  largely  concerned  with  the  bundle 
(T*Q)^ ;  the  fact  that  the  reduced  space  (T*Q)^  still  has  the  factor  seems  to  be  the  reason  the 

connection  on  the  G^  bundle  J“'(p)  (T*Q)^  is  not  sufficient  to  completely  isolate  the 

vibrational  modes  from  the  rotational  ones.  We  believe  that  the  bundle  fills  this  gap.  These 
remarks  aside,  we  turn  to  the  explicit  construction  of  .  To  do  this,  we  first  describe  . 
Recall  that  the  augmented  potential  is  given  by 

V^  =  V  +  L^  (la) 


where 


(lb) 

For  mechanical  systems  undergoing  stationary  rotations  about  i.e.,  G  =  S0(3)  and  G^^  = 
rotations  about  the  axis  which  is  parallel  to  ^),  we  note  that  gives  the  potential  of  the 
centrifugal  force.  Now  define  Ujj^  as  the  subspace  on  which  or  equivalently  looks 
objective  in  the  sense  of  nonlinear  elasticity  (cf.  Marsden  and  Hughes  [1984]).  More  precisely: 
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5.1  Definition 

UiNT  =  {Sqe  T^QK5q,(L,^dl^Xq)>  =  0  forall  ne  and  <5q.X(^cg>=0  foran  v)  ® 

where  the  first  pairing  is  the  natural  pairing  between  vectors  and  one  forms  while  the  second  is  the 
metric  inner  product. 

Since  has  a  critical  point  at  (by  the  principal  of  symmetric  criticality)  and  V  is  G- 
invariant,  we  find  that 

<5q,  (L^^dC^)(q)>  =  S^V^fq^XSq.  TigCq))  (3) 

and  so  we  see  that  the  geometric  condition  (2)  is  exactly  what  is  needed  to  block  diagonalize 
5^V^(q^)  within  S.  In  coordinates,  the  first  condition  on  5q'  defining  Ujjgj  is  the  geometric 

condition 


Sq'  T  e  ^  [Ai  ^  (A,'  A-  g,„)]  =  0;  (2-) 

the  second  condition  is  just  the  defining  condition  on  S.  Now  we  are  ready  to  define  . 

5.2  Definition  =  {5z€  T^T*Q|6q€  Ujjjt  and  5ze  ker[TJ(z^)]}  cj.  (4) 

5.3  Proposition  S  =  Jrig  ®  ^int 

This  is  easy  to  check.  The  idea  is  that  Jrjq  n  =  (0) ,  that  dim  =  dim(g/g^)  and 

that  is  determined  by  dim(g/g^)  equations.  Also,  we  write 

'%rr  -  '^iNT  ® 

where  =  {5p  -  A^q^)  j  5p  6  U^j  is  the  dual  space  with  a  momentum  shift  by  A^  (see 
equation  (10)  of  §4).  The  relation  (5)  is  really  a  coordinate  description;  to  do  it  intrinsically,  we  use 
the  metric  connection  to  split  T^T^Q  =  T^^  Q  ©  (T^^  Q)*  (this  split  is  in  fact  nothing  more  than 

what  we  do  in  coordinates  to  identify  accelerations  and  momenta  with  vectors)  then  we  take  the 
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horizontal  and  vertical  splitting  of  Ujj^  in  with  the  vertical  component  followed  by  the 
momentum  shift  by  A^(q^). 

We  remark  here  that  even  if  G  is  abelian  (for  instance,  G  =  in  the  case  of  planar 
coupled  rigid  bodies)  then  the  decompositions  are  not  triMal;  while  =  {0}  in  this  case, 

=  is  still  not  a  trivial  decomposition. 

Now  +  (Pg,  and  we  have  arranged  for  to  be  block  diagonal.  As  far 

as  is  concerned,  we  compute  in  coordinates  that 

i  g‘J(P|  -  A”4'’) .  (6) 

Thus,  since  pj  =  gjj^  A^  at  equilibrium,  we  get 

5^K^(Zg)  •  (6z,  5z)  =  g'J  8pj  5pj  .  (7) 

It  is  clear  that  block  diagonalizes  from  within  by  construction. 

Regarding  the  block  diagonalization  of  on  5rjq  ©  ,  we  shall  use  some  further 

interesting  identities. 

First,  here  is  an  equivalent  characterization  of  Ujjgj  in  terms  of  superposed  motions: 

5.4  Proposition  Let  q^  €  Q  be  a  curve  tangent  to  5q  ar  q^  ,  let  q  6  and  let  = 
Ade,p(£^)q  .  Then  Ujjyu  is  characterized  by  those  5q  orthogonal  to  Tq^(G^^  ■  q^)  and  satisfying 

^<^(q£).me)Q(qe)>|^_^  =  0  (8a) 

or,  equivalently, 

^<^Q(exp(e4)qg) ,  qQ(exp(e4)q^))  =  0.  (8b) 

This  is  verified  by  a  direct  coordinate  calculation.  We  can  lift  this  expression  to  get  an 
alternative  characterization  of  .  We  consider  the  momentum  map  J  restricted  to  g^^  and 

regarded  as  a  function  on  TQ.  In  other  words,  for  ^  €  g^^  ,  set 


J(C)(5q)  =  <CQ(q).5q)  =  A- C*(5q)j. 


(9) 
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Now  consider  the  condition 


|j(0(6q)|^_^=  0  (10) 

where  ^  is  to  evolve  as  ^  =  [4,  CJ  which  is  consistent  with  (8a)  and  ^  e  ;  here  4  is  the  Lie 

algebra  element  giving  the  relative  equilibrium.  Equation  (10)  defines  a  condition  on  T(TQ).  We 
shall  regard  it  as  a  condition  on  T^(T*Q)  via  the  Legendre  transform.  For  simplicity  we  still 

write  the  resulting  condition  as  J  =  0. 

5.5  Proposition 


~  fi)  ^  S- 


(11) 


The  condition  TJ(Zg)  •  5z  =  0  reads 


3A* 

5Pi  A^(q)  +  pj  ^  dq*'  =  0 


(12) 


and  using  this,  one  can  express  the  conditions  defining  entirely  in  terms  of  5q.  This 
recovers  the  space  Ujj^ ,  which  models  T|qj(Q/G),  and  then  one  gets,  as  before. 


'%rr  ~ 


Uint®  ^4^ 


§  6  Block  Diagonalization 

The  block  diagonalization  results  for  follow  from  two  basic  formulas: 
6.1  Proposition  Let  Az  e  5riq  and  5z  e  T^  P.  Then 

6^H^(Zg)(Az,  5z)  =  ^(CQ(q).  5q>- ([^,  q],  5j(z^)-5z) 


(1) 
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where  Az  has  associated  "n  and  C  as  in  (3b)  and  (7)  of  §4. 

6.2  Proposition  Let  6zj  and  e  ;  then 

8^H^(Ze)(5zi,  5z2)  =  S^K^CZg)  •  (5zj,  6x2)  +  62v^(q^)(5qj,  Sqj)  (2) 

Proposition  6.1,  which  is  proved  by  direct  calculation,  shows  that  8^H^(Zg)  block 
diagonalizes  on  ,  i.e.,  if  Az  €  and  5z  e  then 

52h^(z^)(Az,  5z)  =  0  .  (3) 

Proposition  6.2  then  follows  from  our  earlier  calculations.  It  also  follows  that  if  Az  e  5jyG  and 
A  z  €  then 


52h^(z^)(Az,  Az)  =  ^(CgCq).  Cq>  (4) 

which  is  a  generalization  of  the  rigid  body  second  variation  formula  for  motion  on  the  coadjoint 
orbit  with  the  metric  We  summarize: 

6.3  Theorem  The  relative  equilibrium  is  formally  stable  (with  6^H^(Zg)  on  S  positive 

definite)  if  and  only  if 

*  CgCq))  S^jc 

and  ii  8^V^(q^)  is  positive  definite  on  Ujfjj 

We  note  that  S^V^Cq^)  separates  (in  coordinates  on  Ujjjj)  into  5^V(qg)  plus  a  term 

quadratic  in  Thus,  ii  is  implied  by  a  condition  of  the  form  ||  4  II  ^  V^inin,  where  ||  ||  is  a 
suitable  norm  and  is  the  minimum  (non-zero)  eigenvalue  of  5^V(q^);  one  has  to  take  care 
here  that  V  itself  does  not  have  a  critical  point  at  q^,  so  5^V(q^)  does  not  make  intrinsic  sense. 
To  see  how  this  works  in  examples,  see  Simo,  Posbergh  and  Marsden  [1989]  and  Lewis  and  Simo 
[1989]. 


As  far  as  the  symplectic  form  is  concerned,  we  have 
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6.4  Theorem  Let  Az  €  5^0  and  6z  €  P.  Then 

Q(Zg)(Az,  5z)  =  -  <11, 6J(z^)  •  8z>  +  (CgCOe),  8q>.  (5) 

Notice  that  in  an  appropriate  sense,  8^H^(Zg)  on  ^rjgXT^P  is  the  time  derivative  of  the 
symplectic  fom  ft! 

From  (5)  and  (7)  of  §4  one  finds  that  on  5rjg  x  Jriq,  ft  gives  the  coadjoint  orbit 
symplectic  form 


ft(z^)(Az,Az)  =  -<^^,^1^,n]>. 


(6) 


while  on  x  we  have  the  cross  terms 


ft(z^)(Az,  8z)  =  (Cgfac).  5q>, 


(7) 


which  depend  on  the  8q  components  alone. 

We  can  summarize  the  situation  with  the  following  matrices: 


s/V 


^INT 


52H^(z^)  = 


Generalized 
Rigid  Body 
Second  Variation 


0 

0 


for  the  amended  energy,  and 
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^INT 


Coadjoint  Orbit 

Internal- 

Symplectic  Form 

Rigid  Coupling  (7) 

Q  = 


Internal- 

Rigid  Coupling  (7) 


Canonical  symplectic  form 
plus  a  "magnetic"  term 


0 


For  information  on  the  "magnetic"  term,  and  its  interpretation  as  a  curvature,  we  refer  the  reader  to 
Kummer  [1981].  Also,  the  coupling  terms  can  be  interpreted  in  terms  of  the  curvature  of  the 
connection  on  the  coadjoint  bundle  T*Q  T*(Q/G)  ;  see  Montgomery  [1986]  and  Lewis, 
Marsden,  Montgomery  and  Ratiu  [1986]. 
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